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Summary 

The viscous flow Field near the surface of a hovering 
rotor blade was studied for blade twist distributions typical 
of a tilt rotor blade and a conventional helicopter rotor 
blade. Three blade geometries were studied, including a 
tilt rotor blade twist distribution (baseline), conventional 
helicopter rotor blade twist distribution, and the baseline 
twist distribution with 2 deg of precone. The results give 
insight into the delayed stall phenomenon often observed 
for highly twisted rotors. Calculations were performed for 
a high thrust condition near stall using the thin-layer 
Navier-Stokes CFD code TURNS. Effects of built-in 
twist on section force coefficients, skin friction, 
velocities, surface pressures, and boundary layer shape 
factor are discussed. Although the rotor thrust coefficient 
was nominally the same for the cases using the two twist 
distributions, large differences were found in the section 
in-plane and normal force coefficients. These preliminary 
results imply that the blade outboard region, rather than 
the inboard region, provides the majority of the 
performance advantage of the baseline case over the low 
twist case. Skin friction, velocities near the blade, and 
surface pressures for the two twist distributions reveal 
significant differences in the blade outboard region. 


Notation 

aoo freestream sound speed 

c chord 

Cj-ef reference chord, chord at 0.75 radial station 

(=25.075 in) 
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Cn section normal force coefficient, 

section normal force/CO.Spoofftr)^) 


Cp surface pressure coefficient, 

(p-poo)/(0.5poo(tlr) 2 ) 

Ct section thrust coefficient, 

section thrust/(0.5poo(flR)^c re f) 

Cx section in-plane force coefficient, 

section in-plane force/(0.5p 0 o(£2R) 2 c r ef) 

CT rotor thrust coefficient, 

rotor thrust/(7cR^p(QR) 2 ) 

j,k,& grid indices 

tc grid point in blade normal direction defining 
boundary layer edge 

H boundary layer shape factor, 8*/0 

M t j p blade tip Mach number, £2R/aoo 

Poo free stream static pressure 

r blade radial coordinate, radial location/R 

R blade radius (=19 ft) 

Re Reynolds number, £2Rc re f/v 

x,y,z,t inertial coordinates 

u,v,w nondimensional inertial space velocities 

ub» v B> w B u,v,w rotated from (x,y,z) directions to 
(£ji,Q directions 



8 boundary layer thickness 

8* boundary layer displacement thickness 

v kinematic viscosity at sea level, 1.5723 xl(H 
ft 2 /s 

p density 

poo free stream density 

0 momentum thickness 

Q rotor speed (381 rpm) 

blade-fixed computational coordinates 


Introduction 

Interest in the flow field of rotors in hover or axial 
flight has a long history. In particular, the demonstrated 
ability of highly-twisted rotors or propellers to delay stall 
and continue to generate lift under high thrust conditions 
has prompted numerous studies over the years. The 
ability of rotating, high incidence blade sections to 
maintain attached flow when non-rotating sections under 
similar conditions undergo flow separation has generated 
several explanations for this phenomenon, such as 
outward flow of the boundary layer, Coriolis effects, and 
rotor wake contraction. 

One of the earliest experiments addressing the delayed- 
stall problem is Himmelskamp (Ref. 1). The 
Himmelskamp experiment investigated a two-bladed 
airscrew in axial flow in a duct. High section lift 
coefficients were attributed to the thinning of the boundary 
layer caused by the radial flow produced by the centrifugal 
force. Himmelskamp noted that the Coriolis force also 
contributes to keeping the flow attached at high angles of 
incidence. After the work of Ref. 1, several theoretical 
studies of laminar boundary layers on rotating flat plates 
emerged (Refs. 2 - 5). Harris (Ref. 6) summarized these 
theoretical efforts and also reviewed the experimental work 
completed through the mid- 1960’s in his study of radial 
flow effects on rotor blades. 

At the time of Ref. 6, little or no work had been done 
on turbulent boundary layers on rotating blades. Tanner 
and Yaggy (Ref. 7) conducted an experimental 
investigation, using a chemical film technique, of the 
boundary layer on a UH-1 untwisted tail rotor blade in 
hover and compared results with non-rotating wind tunnel 
measurements. Reference 7 did not find any evidence of 
outward or radial flow of the boundary layer for the 
rotating case. Also, no apparent change in the boundary 


layer thickness was found when comparing rotating and 
nonrotating cases. However, Ref. 7 was not able to 
explore stalled conditions and suggested that outward flow 
of the boundary layer was still possible for regions of 
separated flow. McCroskey and Yaggy (Ref. 8) and 
Dwyer and McCroskey (Ref. 9) provided analyses for a 
laminar boundary layer on a rotating flat plate in forward 
flight, which led to the analysis of a turbulent boundary 
layer on a rotating flat plate by McCroskey, Nash, and 
Hicks (Ref. 10). 

By the 1980’s, literature from the wind turbine 
industry also began appearing on the subject of delayed 
stall, for example, Milborrow (Ref. 11) and Wood (Ref. 
12). Also during this time period, hover tests of two full- 
scale tilt rotors (Refs. 13-14) and a 0.658-scale V-22 rotor 
(Ref. 15) were conducted. Narramore and Vermeland (Ref. 
16) made one of the first attempts to address the delayed 
stall phenomenon using computational fluid dynamics 
(CFD) methods. Reference 16 performed calculations for 
the 0.658-scale rotor of Ref. 15, showing that the flow 
remained attached at inboard blade regions for very large 
collective angles while outboard regions experienced 
separated flow. No evidence of radial flow was found. 
However, the CFD analysis relied on a separate lifting- 
surface code to provide the wake induced angle-of-attack 
distribution at the blade. In addition, the grid spacing used 
near the blade surface was rather coarse. Tung and Branum 
(Ref. 17) conducted a hover test of a model tilt rotor with 
surface pressure measurements and compared these data 
with hover free wake codes coupled with a full-potential 
rotor code. The codes were unable to predict the pressure 
coefficient distributions for radial stations inboard of 50 
percent. Tsung and Sankar (Ref. 18) investigated the 
effects of rotation on flow separation using CFD. 
However, the blade modeled was an untwisted blade with a 
BERP tip and results were only provided for the outer 35 
percent of the blade. Calculations revealed that blade 
rotation allowed the flow to remain attached at incidence 
angles where separation occurred on a fixed wing. Most 
recently, Felker (Ref. 19) used five analytical models to 
predict tilt rotor hover performance and compared these 
predictions with the data from Refs. 13-15. In general, 
predicted power was greater than the measured power, 
especially at high thrusts. Several of the analyses relied 
on two-dimensional airfoil data for lift and drag, resulting 
in premature lift loss when the calculated angle of 
incidence at the blade section was beyond the two- 
dimensional stall angle. 

Despite many studies investigating the delayed stall 
phenomenon, a complete and consistent explanation for 
this difficult problem does not exist. Accurately 
calculating a three-dimensional, rotating, turbulent 
boundary layer with areas of flow separation may not be 
possible for some time. Experimental measurements of 
the boundary layer under these conditions is also an 
intimidating task. The current CFD technology, however, 
does provide a means for scrutinizing the flow field details 
and may provide insight into the stall delay phenomenon. 
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The objective of this study is to shed light on the 
delayed stall phenomenon by examining the viscous flow 
field near the surface of a hovering rotor blade. 
Calculations are performed using a thin-layer Navier- 
Stokes CFD code. The paper will discuss the flow solver 
and solution, the grid geometry, the method of computing 
the boundary layer parameters, and the assumptions and 
approximations used in this study. The effect of blade 
twist and precone on the flow are presented. Two twist 
distributions, representative of a tilt rotor and of a 
conventional helicopter rotor blade, are used. Although 
precone has not been modeled in previous studies, precone 
does affect the velocities and centrifugal force seen by the 
blade. Thus, the precone effect is studied here. By 
examining the section force coefficients, regions of the 
blade are identified which warrant closer study. Then, the 
effect of blade twist on boundary layer parameters such as 
skin friction vectors and shape factor, in addition to 
surface pressure distributions are discussed. Velocity 
profiles for different chordwise and spanwise locations are 
also presented. 


Description of Flow Solver and Solution 

TURNS (Transonic Unsteady Rotor Navier-Stokes) is 
a versatile and robust Euler and Navier-Stokes CFD code 
developed by Srinivasan, et al (Refs. 20-23) at NASA 
Ames Research Center under U. S. Army Research Office 
support. This code has been used in many investigations 
of rotorcraft aerodynamic and acoustic phenomena, 
examples of which can be found in Refs. 20-23. 
Srinivasan, et al (Ref. 21) provides details of the finite 
difference scheme. A brief description of the governing 
equations and numerical scheme, taken from Ref. 23, is 
provided below. The governing equations are given by 


dtQ + + a^F +a^G = -L a^s (i) 

Re 

where 

Q = [ p pu pv pw e ] T , Q = Q/J (2) 


E, F, G are the convective flux vectors and S is the 
viscous flux vector using the thin-layer approximation. 
The blade-fixed computational coordinates are defined by 
(£,T|,Q. The flux vectors have been normalized by the 
Jacobian J. The density, mass fluxes, and energy per unit 
volume are nondimensionalized by the far-field reference 
quantities. Re is the Reynolds number based on Mtip and 
a reference chord (defined as the chord at 0.75R for this 
investigation). The reference chord and sound speed are 


the characteristic length and velocity scales, respectively. 
The equation of state for a perfect gas and Eq. 1 describe 
the entire flow field. The numerical method uses a Roe 
upwind-biased scheme for all the coordinate directions 
with higher-order MUSCL-type limiting for the right-hand 
side. A lower-upper Gauss-Seidel implicit operator is 
used for the left-hand side. The governing equations are 
solved in the blade-fixed coordinate system, requiring a 
transformation that introduces source terms that account 
for the centrifugal force produced by blade rotation. The 
space metrics are evaluated using a finite volume 
formulation, while a finite difference formulation is used 
for the time metrics. The turbulence model in TURNS is 
the Baldwin-Lomax model. 

The boundary conditions described in Refs. 22-23 are 
repeated here. At the blade surface, a no-slip condition is 
used with a finite velocity due to blade rotation. The 
surface pressure is obtained through the normal 
momentum equation and the density is evaluated by the 
adiabatic wall condition. In the azimuthal direction, a 
periodic flow condition is used, hence the need to only 
model one blade. At the blade root plane, the inertial 
velocity components are obtained by assuming the 
velocity gradient in the spanwise direction is zero at the 
root. This boundary condition can affect the development 
of the root vortex and requires improvement. To prevent 
flow recirculation in the computational domain, the far- 
field flow is determined by superimposing the flow from a 
three-dimensional sink (placed at the center of rotation) 
and momentum theory. Srinivasan, et al (Ref. 23) 
provides expressions for the velocity entering the 
computational domain at the upper and side boundaries and 
the exit velocity at the lower boundary. Since total 
performance predictions are not the objective of this study, 
the far- field boundary conditions are considered adequate. 

The solution is provided by TURNS in terms of Q at 
every grid point. For this study, Q, the (x,y,z) grid, and 
the metrics at each grid point calculated by TURNS are 
then used by a separate program to compute the skin 
friction and inertial velocities (rotated from the (x,y,z) 
directions to the (^,tj,Q directions), the boundary layer 
shape factor, and surface pressure coefficients. 


Description of Grid 

The finite-difference grid is of C-H topology, with C 
in the wrap-around direction and H in the spanwise 
direction. The grid dimensions are 181 points in the 
wrap-around direction, with 49 points in the spanwise and 
normal directions. The outer boundary is one rotor 
diameter above and below the rotor plane. Only one blade 
is modeled. The distance to the first point off the blade 
surface in the normal direction is 0.00004c re f which is 
approximately 0.27% of the boundary layer thickness at 
the 0.75 radial station. The grid stretching in the normal 
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direction gives 15-20 grid points in the boundary layer. 
The spanwise location of the outer boundary is one rotor 
radius beyond the blade tip. Figure la and lb show the 
computational coordinate system and a plan view of the 
grid geometry, resepcdvely. The inner boundary is at the 
blade root (0.23 radial station) and the reference chord is 
located at the 0.75 radial station. 

For this investigation, the grid generator was 
modified to incorporate twist distributions typical of tilt 
rotor blades and conventional rotor blades. In addition, the 
option of blade precone was included. The airfoils used 
are similar to those developed at Bell Helicopter Textron. 
Inc.; all blade geometries incorporated the same spanwise 
distribution of airfoils. Figure 2a and 2b show the two 
twist distributions studied in this paper, without and with 
the collective pitch, respectively. The collective pitch is 
16 deg for the high twist case and 19 deg for the low twist 
case. Different collectives were required in order to match 
rotor thrust coefficients. The high twist is representative 
of a tilt rotor blade and is defined as the baseline case for 
this study. Figure 3 shows spanwise cuts at the blade 
root for die two twist distributions. The collective pitch 
is included in the geometry. In Fig. 3a, the severe twist 
and a requirement by the grid generator to enforce a zero- 
slope at the outer boundary cause a kink in the grid aft of 
the trailing edge. The grid is designed to capture the flow 
at the blade surface, but is not optimized to capture tip 
vortex formation, the trailed wake behind the blade, or the 
wake underneath the blade from the preceding blades. 
Since performance is not the focus of this investigation, 
the grid is considered sufficient for analyzing the (low on 
and very near the blade surface. Applications such as Ref. 
23 imply that the boundary conditions and grid are 
adequate to define the aerodynamic environment near the 
blade. 


Definition of Rotation Matrices 

In order to calculate skin friction vectors and shape 
factor values on the blade surface, the inertial space 
velocities and the stress tensor must be rotated to the 
directions. Orthogonal metrics x^, £ x , etc., can be 
used to form rotation matrices for this transformation. 
However, because the grid is only nominally orthogonal 
at the blade surface for the geometries in this study, the 
rotation matrices defined by the metrics will not be 
orthogonal matrices. Therefore, rotation matrices are 
developed directly from the surface geometry. This is 
done by first constructing a vector tangent to the £- 
direction at the blade surface for a (j,k) point, then making 
t£ into a unit vector, u^. Next, a vector e rj tangent to the 
q-direction is used to construct a vector orthogonal to 
u^. The vector tq is then made into a unit vector u^. 
The third orthogonal vector, u^, is obtained by taking the 
cross product between and urj. Finally, the rotation 


matrices Cba and Cab (between inertial space A and the 
blade surface B) are constructed. The procedure is shown 
in the following equations. 


0-5 (*j+l,k,l * *j-l,k,l) 




0.5 (yj+i,k,i - yj-i.k.i) 


Ctl = 


0.5 (Zj + i,k.l - z j-l,k,l) 

0.5 (Xj.k+l,l - Xj.k-l.l) 
0.5 (yj.k+i,i - yj.k-i,i) 
0 5 (zj.k+ 1,1 - z j.k-l.l) 


*ti = <[i] - U ^ T ) £ -n 


(3) 


(4) 

(5) 



It 1 

-i- , X Utl 


( 6 ) 


•V 


[Cba] = 


[Cab]=[ “il “C ] 


(7) 


( 8 ) 


The rotation matrices are defined at the blade surface (6=1) 
for a particular (j,k) location and then used for all 6 at that 
location. The velocities and stress tensor are then 
transformed by 


UB 

= [Cba] 

U 

VB 

V 

. wb _ 


. w . 


■ 

x kr\ 

U-T 

= [Cba] 

^xx 

T *y 

^xz 

x ll^ 

Ititi 


*yx 

T yy 

'tyz 


T Sn 



^zx 

^zy 

^zz 


[Cab] 

( 10 ) 


Although TURNS assumes a thin viscous layer near the 
blade surface to compute Q, the complete expressions for 
the tensor elements in the (x,y,z) system (Eq. 10) were 
evaluated before being transformed. 
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Determining Boundary Layer Parameters 

The following equations from White (Ref. 24) were 
used to calculate the displacement thickness, momentum 
thickness, and shape factor of the boundary layer at each 
(j Jc) location on the blade upper surface. 


The skin friction values in the and q-directions are 


computed using the following: 


Cf % = T^/(0.5poo(Qr) 2 ) 

(17) 

Cf n<; = Xt,j;/(0.5poo(£2r) 2 ) 

(18) 


8*= f o--e-jL)dz= £ (i-P£- >Ae (id 
Jo P eU e (T, % Ufa 


£._u_(i.jL)dz= y (PL 
0 Pe U e Ue e=l Pk 


ue* ueo 


(12) 

H = £l (13) 

e 

Ae = t ( x e+i - x e ) 2 + (ye+i - ye ) 2 + ( z e+i - z e ) 2 i 05 ( 14 ) 

Uec = [(“Bee - u Bl) 2 + ( v Bee - VBl) 2 1°* ( 15 > 


Since a definition does not exist for the boundary layer 
edge of a rotating twisted blade, an approximation using a 
flat plate analogy is used. The edge of the boundary layer 
was first estimated by calculating the boundary layer 
thickness, 8, for turbulent flow over a fiat plate using the 
Reynolds number at the 0.75 radial station. Next, the 6- 
value (£max) corresponding to a distance 8 away from the 
blade surface was determined for the 0.75 radial station and 
used as an upper limit on te (the t -value defining the 
boundary layer edge). The boundary layer thickness will 
vary for different (j,k) locations on the surface; therefore, 
It was determined by searching for the maximum velocity 
magnitude from t=l to £^max at each (j,k) location. 
This maximum velocity was then defined as U£e a °d the 
corresponding € -value defined as 6e. The dot product of the 
edge velocity vector and the velocity vectors from C=1 to 
€=Ce was then computed to obtain the components u£ in 
the direction of the edge velocity, as shown in the 
equation below (note that the velocity, u£, is zero at t=l): 


(u Be : u Bl) (UB £e - u Bl) + (vBe - v Bl) ( v BCe 1 ^l) 

UU 


06 ) 


Finally, 8*, 0, and H were computed. 


Results 

The effects of blade twist and precone on section force 
coefficients and flow parameters are studied by performing 
calculations for three blade geometries: tilt rotor blade 
twist distribution (baseline case), 10 deg linear twist 
distribution (low twist case, simulating a helicopter case), 
and tilt rotor blade twist distribution with 2 deg precone. 
Figure 2 shows the baseline and low twist distributions 
used for the calculations. All three geometries incorporate 
the same spanwise airfoil distribution and chord 
distribution (chord distribution shown in Fig. 4). For the 
three cases, Mtip = 0.68, Reynolds number based on Mtip 
and the reference chord =10.1 million, and CT=0.017 
(nominally). The collective pitch for the baseline twist 
cases is 16 deg; the low twist case collective is 19 deg 
(Fig. 2b). Tabic 1 provides a summary of the rotor 
parameters and operating conditions for the calculations 
presented. 

Residual and Rotor Thrust Coefficient 
Behavior 

TURNS was run for 7000 iterations for each of the 
three cases. Figure 5 shows the residual behavior. The 
residual for iterations less than 1000 were computed using 
a grid coarsened by using every other grid point in all 
three directions. In Fig. 5, RSUM (L2 norm) represents 
the residual over the entire grid volume. Less than two 
orders drop in magnitude of the L2 norm is shown in the 
first 2000 iterations for all three cases. Since the largest 
contributor to the L2 norm was usually located at the 
outer boundary, a residual within the grid volume of 
interest was computed. RSUMX represents the residual 
over a volume defined by the blade surface from the root 
to 0.96 radial station and out to 30 grid points away from 
the surface (a distance of approximately 0.10-0.15 of the 
reference chord). RSUMX shows several orders of 
magnitude drop in the first 2000 iterations. 

The rotor CT behavior is shown in Fig. 6. The 
oscillations are possibly caused by the movement of the 
near field tip vortex trajectory. These oscillations are 
more severe if the grid near the tip becomes coarser. The 
CT value averaged from 2000 to 7000 iterations is 
nominally 0.017 for all three cases. 


3.3- 5 



Sectional Force Coefficients 

Examination of the spanwise behavior of the 
sectional forces is instructive for identifying flow regions 
of interest on the blade. Figure 7a and 7b present the 
section thrust (Ct) and in-plane (Cx) force coefficients, 
respectively. Figure 7a shows that the overall thrust 
produced by all three cases is nominally the same. The 
baseline twist case (with and without precone), however, 
produces more thrust for radial stations less than 0.6 
compared to the low twist case. Outboard of 0.85 radial 
station, the thrust is nearly the same for both twist cases, 
even though the low twist case has a larger blade pitch in 
this region (Fig. 2b). Figure 7b shows that the Cx 
distribution is markedly different between the two twist 
cases, especially outboard of the 0.6 radial station. Since 
Cx times the blade radial station integrates to the rotor 
torque, the differences in Cx between the two twist cases 
at the outboard region is more important than the 
differences at the inboard region when computing rotor 
performance. The higher Cx near the blade tip for the low 
twist case implies a performance worse than the baseline 
twist case, which is consistent with previous tilt rotor 
measurements. The low twist case has a higher geometric 
pitch than the baseline case outboard of 0.6 radial station 
which will contribute to a higher induced drag. However, 
the magnitude of the wake effect compared lo the increase 
in profile drag is unknown. 

The behavior of the blade normal (Cn) force 
coefficient based on the local chord and radial station is 
also revealing (Fig. 8). Cn is normal to the local blade 
chord. The baseline twist blade sections generate more 
normal force than the low twist blade for radial stations 
less than 0.6. The Cn begins decreasing at 0.8 radial 
station for the low twist blade; this drop-off in Cn is 
delayed until approximately 0.85 for the baseline twist 
blades. If Cn is assumed to be approximately the lift 
coefficient, some comparisons can be made with two- 
dimensional measurements. For example, the maximum 
2D lift coefficient for the airfoils used between 0.2 and 
0.45 radial stations is 1.2- 1.3. Clearly, Fig. 8 shows the 
baseline twist cases exceed this value, while the low twist 
case does not This behavior is similar to that seen in 
other studies in the literature (Refs. 17 and 19). From 0.4 
to 0.85 radial station, however, both twist distributions 
give Cn values below the maximum 2D lift coefficient 

Since the above results indicate the blade outboard 
region has a larger impact on performance than the inboard 
region, as indicated by the dramatic differences in Cx in 
Fig. 7b, a closer look at the blade region outboard of 0.6 
radial station is warranted. 

Skin Friction 

Figures 9-11 show skin friction vectors on the blade 
upper surface outboard of 0.789 radial station for the 
baseline twist low twist, and baseline twist with precone, 
respectively. For clarity, every other chordwise grid point 


is shown. For radial stations inboard of 0.789, the skin 
friction vectors are essentially chordwise for all three cases 
without any flow separation. There is a dramatic 
difference between the baseline and low twist cases 
(compare Figs. 9 and 10). The low twist case shows a 
significant amount of outward spanwise flow, starting 
from the 0.789 radial station. The baseline twist case 
shows the flow is chordwise out to 0.836 radial station. 
The 0.789 and 0.836 radial locations correspond to the 
approximate spanwise locations where Cn (Fig. 8) begins 
to decrease for the low and baseline twist cases, 
respectively. The addition of 2 deg of precone to the 
baseline twist does not have a significant effect (compare 
Figs. 9 and 11). Both twist cases appear to have a 
counterclockwise flow pattern over a portion of the blade: 
between 0.857 and 0.908 radial stations for the baseline 
twist and between 0.789 and 0.892 radial stations for the 
low twist case. The velocity vectors (with components 
ub-ub 6=1 Bnd v B6 v B6=l) at approximately 0.01 chord 
distance above the blade are shown in Figs. 12 and 13 for 
the baseline and low twist cases, respectively. The 
velocity vectors indicate there is some spanwise flow near 
the blade surface and also some regions of swirl; the 
velocity vectors are consistent with the skin friction 
vectors. 

Velocity Profiles 

Velocity profiles in the boundary layer are examined 
for several chordwise stations at the 0.857 radial station 
for the baseline and low twist case in Figs. 14 and 15, 
respectively. The chordwise, spanwise, and normal 
velocity components are ub-ub£= 1 , vbC _v b 6=1. and 
wb£"Wb6=i. respectively. In Fig. 14, the flow remains 

attached over the blade chord. There is a small 
contribution from the spanwise component; the normal 
component is essentially zero. For the low twist case, 
Fig. 15 indicates that the flow is already separated close to 
the leading edge, after which the chordwise velocity 
becomes increasingly negative (reverse flow). The 
spanwise velocity grows in magnitude as the trailing edge 
is approached until the flow is almost entirely spanwise, 
which is consistent with Figs. 10 and 13. 

Pressure Coefficient Profiles 

The surface pressure coefficient (Cp) distributions cor 
responding to the radial station (0.857) of Figs. 14 and 15 
is shown in Fig. 16. The baseline case with precone has 
nearly the same Cp distribution as the baseline case 
without precone. In Fig. 16, the lower surface Cp 
distributions are essentially the same. On the upper 
surface, the low twist case has a severe adverse pressure 
gradient starting from the leading edge and thereafter 
remains relatively constant until the trailing edge. This 
behavior indicates separated flow (seen in Fig. 15) over 
most of the upper surface. The difference in the upper 
surface Cp distribution contributes to the difference in Cx 
at this radial station (Fig. 7b). 
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Shape Factor 

The shape factor is a useful parameter for identifying 
separation points for two-dimensional boundary layers. 
For turbulent flat plate flow, the shape factor is on the 
order of 1.2 (Ref. 24), with separation occurring near H=2 
(Schlicting, (Ref. 25)). Computing shape factors for 
three-dimensional, rotating blades at large incidence angles 
using Eqs. 11-13 may be inappropriate, especially using 
the approximate method for finding the boundary layer 
edge. The intent here, however, is to ascertain whether 
the behavior of H, not necessarily the value of H, can 
serve as another indicator of the flow behavior. Figure 17 
presents chordwise shape factor distributions for several 
spanwise stations for both the baseline and low twist 
cases. Comparing Fig. 17 with the skin friction plots 
reveals the chordwise location where H undergoes a sign 
change gives an indication of the change in direction of 
the skin friction vectors (from outward spanwise to inward 
spanwise). The sign change also reflects the reversal in 
sign of U£ e > the velocity at the boundary layer edge. 


Concluding Remarks 

The objective of this study was to shed light on the 
delayed stall phenomenon by examining the viscous flow 
field near the surface of a hovering rotor blade. The 
Navier-Stokes CFD code TURNS was used to perform 
calculations for three blade geometries which included a 
tilt rotor blade twist distribution (baseline), conventional 
helicopter rotor blade twist distribution, and the baseline 
twist distribution with 2 deg of precone. Calculations 
were performed for thrust conditions near stall. Although 
the rotor thrust coefficient was nominally the same for the 
cases using the two twist distributions, large differences 
were found in the section in-plane and normal force 
coefficients. These results imply that the blade outboard 
region, rather than the inboard region, may provide the 
majority of the performance advantage of the baseline case 
over the low twist case. Skin friction, velocities near the 
blade, and surface pressures for the two twist distributions 
reveal significant differences in the blade outboard region. 

This investigation was entirely computational. 
Given the assumptions of this study, the following 
observations are made: 

1. The in-plane force coefficient is much higher outboard 
of 0.6 radial station for the low twist case compared to the 
baseline twist case. This implies a performance penalty 
for the low twist case. 

2. The normal force coefficient is much higher inboard of 

0.6 radial station for the baseline twist case compared to 
the low twist case. 


3. For radial stations outboard of 0.789, there is more 
outward spanwise flow on the blade upper surface for the 
low twist case than for the baseline twist case. There is 
evidence of swirling flow near the blade surface for both 
twist cases. 

4. For radial stations inboard of 0.789 radial stations, the 
upper surface flow remains attached and chordwise in 
direction for both the baseline and low twist cases. 

5. The sign change in shape factor provides an indication 
of spanwise flow direction. 

6. Adding 2 deg of precone to the baseline twist 
distribution has no significant effect on the blade flow 
field. 

Additional work is needed to elucidate the fluid 
dynamic mechanisms involved in the observed flows. 
Before the above observations can be considered 
conclusive, the following actions are required: improve 
boundary condition at the blade root, smooth kinks in the 
grid caused by high twist, increase grid density, 
investigate the role of the tip vortex and wake on stall 
behavior, investigate alternate turbulence model, and 
determine an equivalent three-dimensional shape factor. 
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Figure 1. Grid geometry, (a) coordinate system, ( b ) planview of 6=1 surface (every other grid point shown) 
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Figure 3. Grid geometry at blade root (every other chordwise grid point shown), (a) baseline twist, 
precone-o deg, (b) low twist, precone -0 deg 
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Figure 4»Spanwise chord variation. 
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Figure 6. CT behavior, a) baseline twist, precone =0 deg , b) 
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Figure 7. S panwise variation of thrust and in-plane force coefficient, (a) thrust coefficient, (b) in-plane 
force coefficient 
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Figure 8. Spanwise variation of normal force coefficients. 
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Figure 9. Skirt friction vectors for baseline twist, precone=0 deg. 
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Figure 11. Skin friction vectors for baseline twist, precone =2 deg. 


3.3 - 20 



































